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The  multiphase  catalytic  decomposition  of  hydrogen  peroxide  into  water  and  oxygen  is  notoriously  sus¬ 
ceptible  to  thermal  runaway  (heat  of  reaction:  -98  kj  mol-1 ).  The  high  surface  area  to  volume  ratio  (S/I/) 
in  a  microscale  packed  bed  (MPB)  reactor  (radius  0.5  mm)  was  investigated  for  reducing  the  risk  of  ther¬ 
mal  runaway  during  hydrogen  peroxide  decomposition  to  oxygen  intended  as  a  fuel  cell  oxidant  aboard 
an  unmanned  undersea  vehicle  (UUV).  A  microscale  reactor  channel  with  a  S/V  of  ~2  x  103  m2  m-3  sim¬ 
ulated  under  convective  cooling  generated  a  significant  heat  rise  (T  rise  ~  100  K),  whereas  a  microreactor 
with  a  higher  S/ V(~ 200  x  103  m2  nrr3 )  achieved  thermal  control  (T rise  <  10  K)  over  the  simulated  reaction 
zone.  Although  thermal  management  was  successfully  accomplished  using  the  higher  S/V,  experimental 
conversions  of  hydrogen  peroxide  to  oxygen  (5-18%)  measured  from  the  outlet  were  lower  than  simu¬ 
lated  conversions  (38-63%).  Simulation  assumptions,  such  as  homogeneously  dispersed  flow  and  perfect 
catalyst  interaction  among  other  factors,  contributed  to  the  discrepancies  between  the  simulated  and 
experimental  degrees  of  peroxide  conversion  to  oxygen.  Even  though  thermal  control  of  the  MPB  was 
achieved,  this  work  indicates  that  mass  transfer  limitations  are  a  factor  in  the  MPB  reactor  during  a  mul¬ 
tiphase  reaction,  like  decomposition  of  hydrogen  peroxide  to  oxygen  and  water,  and  suggests  means  to 
overcome  them  even  on  the  microscale  level. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  major  challenge  facing  the  design  of  air-independent  power 
systems  is  energy  storage.  High  energy  density  and  potentially  low 
operating  costs  make  fuel  cells  an  attractive  option  for  unmanned 
undersea  vehicle  (UUV)  applications  [1,2].  Modeling  of  an  air- 
independent  solid  oxide  fuel  cell  (SOFC)  power  system  for  UUVs 
by  Burke  and  Carreiro  showed  that  an  efficiency  of  45%  is  theo¬ 
retically  possible  based  on  a  2.5  kW  net  output  over  40  h  [3].  To 
compete  with  current  battery  technologies,  fuel  cell  systems  like 
the  UUV  SOFC  must  also  exceed  300  Wh  L-1 ,  underscoring  an  inher¬ 
ent  volume  limitation  associated  with  UUV  power  generation  [4]. 
Additionally,  air-independent  fuel  cell  power  systems  require  an 
oxidant  supply.  The  work  presented  herein  focused  on  supplying 
oxygen  via  hydrogen  peroxide  (H202)  decomposition  and  evaluat¬ 
ing  the  capabilities  and  limitations  of  such  a  microreactor  used  to 
provide  oxygen.  It  is  estimated  that  for  a  2.5  kW  net  power  output, 
a  minimum  total  flow  rate  of  0.6  02  mol  min-1  (14.7  02  Lmin-1 )  is 
required  [3]. 
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H202  readily  decomposes  into  water  (H20)  and  oxygen  (02)  via 
the  reaction  H202(1)  -*  H20(1)  +  l/202(g)  upon  contact  with  a  num¬ 
ber  of  commercially  available  metal  catalysts  [5]  offering  a  dense 
source  of  02  per  unit  volume.  The  availability  of  high  02  density  in 
a  60%  (w/w)  H202  solution  facilitates  ease  of  handling  and  makes 
concentrated  H202  a  valuable  commodity  as  an  air-independent 
fuel  cell  oxidant.  Furthermore,  the  existing  commercial  infrastruc¬ 
ture  of  H202  has  long  added  to  its  appeal  as  an  oxidant  source  [6]. 
Kinetic  experiments  have  shown  that  manganese  dioxide  (Mn02)  is 
a  relatively  inexpensive,  yet  effective  H202  decomposition  catalyst 
with  reported  activation  energies  ranging  from  20  to  44  kj  mol-1 
and  was  selected  as  the  primary  catalyst  for  this  work  [5,7-10]. 
Despite  these  advantages,  it  is  well  established  that  catalytic  H202 
decomposition  is  highly  exothermic  (A H°xn  =  -98kJmol_1)  [11]. 
This  factor  has  historically  limited  its  application  as  a  power  sys¬ 
tem  oxidant,  especially  within  the  U.S.  Navy  community  [12].  In 
addition,  like  other  exothermic  multiphase  reactions,  the  liquid  and 
gas  product  stream  from  H202  decomposition  can  lead  to  undesir¬ 
able  pressure  increases,  irregular  distribution  of  thermal  capacity 
in  conventionally  sized  reactors,  and  the  possibility  of  thermal  run¬ 
away  [13]. 

Microscale  reactors,  possessing  relevant  fluid  dimension  in  sub¬ 
units  <  1  mm,  have  high  surface  to  volume  ratios  (S/V)  resulting  in 
heat  and  mass  transfer  coefficients  capable  of  inhibiting  thermal 
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runaway.  Conventional  reactors  with  a  relevant  flow  geometry  of 
10  cm,  for  example,  have  S/Vs  of  ~10m2  m-3,  whereas  microscale 
reactors  possess  S/Vs  that  are  at  least  two  orders  of  magnitude 
greater  [14].  The  microscale  packed  bed  (MPB)  reactor  used  in 
this  work  included  an  extended  heat  exchange  surface  and  had  a 
S/V  of  1.88  x  105m2m-3.  The  small  characteristic  dimensions  of 
microreactor  subunits  generally  result  in  low  Reynolds  numbers 
and  subsequently  laminar  flow  [15].  Laminar  flow  minimizes  irreg¬ 
ularities  and  stabilizes  flow  in  the  reaction  zone.  Although  laminar 
flow  is  typical  in  many  homogenously  dispersed  product  streams 
within  microchannels,  slug  flow  can  also  develop  in  multiphase 
flow  scenarios  [16].  In  a  study  of  microchemical  decomposition 
of  H202,  two-phase  flow  evolution  of  02  gas  bubbles  and  liquid 
water  from  platinum-  catalyzed  H202  decomposition  in  80-[xm 
diameter  microchannels  occurred  in  three  different  stages:  liquid 
decomposition,  bubble  coalescence,  and  bubble  slugs  [16].  Despite 
the  complexities  of  the  multiphase  product  stream  generated  from 
H202  decomposition,  in  the  present  investigation  we  aimed  to 
demonstrate  thermal  management  feasibility  and  assess  02  pro¬ 
duction  during  the  catalytic  multiphase  decomposition  of  H202 
using  the  MPB  reactor. 

2.  Methods 

2.1  MPB  geometry 

The  semicircular  MPB  channel  was  0.5  mm  deep  and  5  cm  long. 
The  MPB  channel  was  embedded  in  a  rectangular  stainless  steel 
316  block  (7  cm  x  2  cm  x  0.4  cm)  which  acted  as  the  extended  sur¬ 
face  area.  A  Plexiglas  (PMMA)  cover  (7  cm  x  2  cm  x  0.1  cm)  sealed 
the  top  of  the  channel.  Fig.  1  shows  the  dimensions  of  both  the  sim¬ 
ulated  MPB  geometry  (Fig.  1  A)  and  the  experimental  MPB  geometry 
(Fig.  IB). 

The  experimental  reactor  differed  from  the  simulated  MPB 
geometry  slightly  in  that  the  inlets  and  outlets  were  positioned 
differently.  In  the  experimental  MPB,  both  the  inlet  and  outlet 
were  1  cm  from  the  MPB  edges.  The  inlet  of  the  simulated  reac¬ 
tor  microchannel  was  2  cm  from  the  edge  of  the  MPB  stainless  steel 
block  and  the  outlet  of  the  simulated  reactor  microchannel  was 
flush  with  the  MPB  edge.  Shifting  the  simulated  microchannel,  so 
the  outlet  in  the  model  was  flushed  with  the  external  side  surface, 
eliminated  the  need  for  a  90°  angle  which  would  have  directed  the 
microchannel  towards  the  bottom  external  surface.  This  approach 
minimized  meshing  issues  and  facilitated  model  convergence  nec- 


Fig.  1.  Microscale  packed  bed  (MPB)  reactor  geometry  used  during  simulation  (A) 
and  experimentation  (B). 


essary  for  simulation  of  the  three-dimensional  temperature  and 
concentration  profiles.  The  meshing  algorithm  used  during  sim¬ 
ulation  created  a  mesh  of  over  17,000  tetrahedral  finite  elements, 
which  was  finer  in  the  microchannel  region  where  the  reaction  was 
defined. 


2.2.  Simulated  H202  decomposition 

2.2.1.  Governing  equations  and  assumptions 

A  variety  of  multiphase  flow  regimes,  such  as  homogeneously 
dispersed  gas  in  liquid  or  intermittent  slugs  of  gas  and  liquid,  are 
capable  of  developing  in  microchannels.  The  evolving  flow  depends 
on  a  multitude  of  initial  and  dynamic  flow  conditions.  Based  on  the 
operation  of  the  decomposition  reactor  at  atmospheric  pressure 
and  a  temperature  range  from  approximately  290  K  to  330  K,  it  was 
not  possible  to  determine  a  priori  the  flow  regime  or  regimes  that 
would  evolve  in  the  reactor,  nor  was  the  modeling  algorithm  able 
to  make  this  prediction  [17].  Consequently,  when  solving  the  gov¬ 
erning  equations  for  the  reactor  performing  the  multiphase  FI202 
decomposition,  we  chose  to  describe  the  properties  of  the  materi¬ 
als  in  the  system  as  mass-weighted  averages  of  properties  of  the 
phases  present.  For  example,  the  density  of  the  fluid  within  the 
control  volume  of  the  reactor  was  treated  as  a  single  homogeneous 
phase  whose  density  is  the  weighted  average  of  the  density  of  the 
gas  phase  and  the  density  of  the  liquid  phase  within  the  control 
volume  [18]. 

Mass  balance.  FI202  in  the  vapor  phase  was  considered  negligi¬ 
ble  based  on  the  fact  that  equilibrium  vapor  concentration  of  50% 
(w/w)  IT202  reactant  (the  highest  nominal  concentration  used  in 
this  research)  at  303  K  is  less  than  6%  [19].  Additionally,  the  amount 
of  dissolved  02  in  the  liquid  was  considered  negligible  compared 
to  the  IT202  and  water  liquid  components.  Similarly,  the  amount  of 
water  vapor  in  the  gas  phase  was  considered  negligible  compared 
to  the  02  gas,  provided  the  temperature  remained  below  the  boil¬ 
ing  point  of  water  (373  K).  The  mass  balance  included  the  effects 
of  diffusion  of  FI202  in  a  water-FI202  mixture,  the  bulk  flow  of  the 
evolving  fluid  mixture  (consisting  primarily  of  liquid  FI202  solution 
in  water,  and  02  gas)  down  the  reactor  channel,  and  consumption  of 
IT202  resulting  from  the  decomposition  reaction.  Eq.  (1)  governed 
the  steady-state  mass  balance  of  FI202. 


-DH202eff 


d2cH2o2 

dx2 


~  DH202ett 


92Ch2o2 

dy2 


~  DH202eff 


d2cH2o2 

dz 2 


+  Vx 


9Ch2o2 

dx 


=  rH202 


(1) 


In  Eq.  (1),  DH2o2eff  was  the  effective  diffusivity  of  FI202  in  the 
water-H202  solution.  The  catalytic  reaction  rate  law,  rH2o2,  served 
wholly  as  the  conversion  mechanism.  The  mass,  energy  (Eq.  (2)), 
and  momentum  balances  (Eq.  (3))  were  solved  simultaneously  for 
the  FI202  concentration  (Ch2o2  ),  temperature  (T),  and  velocity  field 
(U)  respectively  using  COMSOL  Multiphysics  3.4.  The  variable  vx 
represented  the  superficial  velocity  of  the  fluid  in  the  x-axis  direc¬ 
tion.  The  rate  of  reaction,  rH2o2,  was  established  assuming  first 
order  kinetics  weighted  via  an  inputted  catalyst  mass  [7].  A  conven¬ 
tional  Arrhenius  relationship  estimated  during  kinetic  experiments 
with  the  Mn02  catalyst  defined  the  reaction  rate  constant  cou¬ 
pling  the  mass  balance  to  the  energy  balance  (Eq.  (2)).  Appendices 
A  and  B  entitled  “Simulation  Constants”  and  “Simulation  Expres¬ 
sions”  specify  all  values  and  expressions  used  in  the  models.  All 
parameters  were  in  SI  units. 

Energy  balance.  The  energy  balance  included  conduction,  advec- 
tive  input  and  output  due  to  the  feed  and  exit  streams,  exothermic 
heat  generation  due  to  reaction,  sensible  heat,  and  the  heat  of  water 
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vaporization.  Eq.  (2)  governed  the  steady-state  energy  balance. 


„  92T  „  d2T  „  92T  9T  Au 

~Keffdx2  -  ^dy2  ~  KeffW  +  Vxf>Cpdx  =  AHrxn%02 

—  Sh~  AHVaPrH2o 


(2) 


fluid  velocity  was  set  equal  to  zero.  The  outlet  boundary  condition 
for  the  momentum  balance  was  defined  by  normal  pressure  and 
flow.  Appendix  C  summarizes  the  simulated  boundary  conditions. 

2.3.  Experimental  H202  decomposition 


Momentum  balance.  The  momentum  balance  was  modeled 
as  a  steady-state,  non-isothermal  flow  using  the  Navier-Stokes 
equations  for  weakly  compressible  flows,  i.e.,  flows  with  Mach 
numbers  <0.3  [20].  Gravity  forces  were  neglected.  Eq.  (3)  gov¬ 
erned  the  momentum  balance,  which  is  a  reduced  form  of  the 
Navier-Stokes  equation  under  steady-state  conditions. 

p(UVU)  =  -Vp  +  r]V2U  (3) 


The  first  term  in  Eq.  (3),  referred  to  as  “convective  acceleration,” 
is  based  on  the  position  derivative  of  the  velocity  field  which  is 
independent  of  time  [21].  The  third  term  describes  the  effect  of 
p,  the  fluid  viscosity,  on  the  velocity  field.  During  laminar  flow  in 
the  x-axis  direction,  the  transverse  components  of  the  velocity  are 
essentially  negligible  and  Eq.  (3)  simplifies  to  Eq.  (3a). 


P 


-Vp  +  r) 


(  d2y*  ,  d2 vx  d2vA 
Y  dx2  dy 2  dz2  J 


(3a) 


Eq.  (4)  is  the  continuity  equation.  Using  low  initial  volumetric  flow 
rates  ranging  from  0.05  ml  min-1  to  0.15  ml  min-1  (corresponding 
initial  fluid  velocities  of  2.12-6.37  mm  s_1)  maintained  the  condi¬ 
tions  defining  weakly  compressible  laminar  fluid  flow. 


V  •  (pU)  =  0  (4) 

The  density,  p,  was  modeled  as  an  average  of  the  liquid  and  gas 
comprising  the  fluid  weighted  by  the  mole  fractions  of  H202,  H20, 
and  02  species  respectively  and  changed  as  the  reaction  progressed 
(Eq.  (5)).  This  average  density  was  applied  under  the  simplifying 
assumption  that  the  generated  02  gas  was  homogenously  dispersed 
throughout  the  fluid  and  the  effluent  stream  maintained  a  well- 
mixed,  dispersed  phase. 


P  =  *h2o2Ph2o2  +  xh2oPh2o  +yo2Po2 


(5) 


2.2.2.  Boundary  conditions 

For  the  inlet  boundary  condition  of  the  mass  balance,  the  H202 
concentration  was  set  equal  to  the  feed  concentration.  The  outlet 
boundary  condition  of  the  mass  flux  was  defined  as  purely  con¬ 
vective,  since  after  leaving  the  reactive  environment  no  diffusive 
effects  occur.  The  microchannel  reactor  walls  constrain  the  total 
normal  mass  flux  to  zero. 

For  the  inlet  boundary  condition  of  the  energy  balance,  the  fluid 
temperature  was  set  equal  to  the  initial  fluid  temperature.  The  flu¬ 
idic  outlet  boundary  condition  of  the  energy  flux  was  defined  as 
purely  convective,  since  the  reactor  ends  at  that  point  and  any 
connecting  tubing  would  have  negligible  conductive  losses  rela¬ 
tive  to  the  main  MPB  reactor  unit.  To  simulate  convective  cooling 
under  ambient  conditions  in  air,  external  heat  flux  defined  the 
boundary  condition  imposed  on  all  external  surfaces  exposed  to  the 
surrounding  environment.  At  the  outer-wall,  air  interface,  which 
occurs  at  the  external  surfaces  exposed  to  the  ambient  environ¬ 
ment,  natural  convection  is  the  limiting  aspect  of  the  dissipation 
of  heat  produced  from  the  reactor  unit.  Continuity  of  both  the  con¬ 
vective  and  conductive  heat  fluxes  defined  the  boundaries  of  the 
internal  surfaces  of  the  microchannel  embedded  in  the  stainless 
steel  block  and  underneath  the  Plexiglas  cover  for  the  energy  bal¬ 
ance. 

In  the  momentum  balance  the  superficial  velocity  was  set  equal 
to  the  initial  superficial  velocity  in  the  x-direction  at  the  inlet.  No 
slip  was  assumed  at  the  MPB  reactor  walls,  and  at  that  boundary  the 


Setup  and  data  collection.  Effluent  02  gas  rates  and  temper¬ 
ature  data  from  the  physical  reactor  were  collected  using  the 
setup  shown  in  Fig.  2.  A  World  Precision  Instruments  SP2001 
syringe  pump  dispensed  nominal  50%  weight  in  water  (measured 
55%,  w/w)  H202  supplied  by  Sigma-Aldrich  at  volumetric  flow 
rates  (VFR)  0.05  ml  min-1  (8.3  x  10-10  m3  s-1)  and  0.15  ml  min-1 
(2.5  x  10-9  m3  s-1 ).  These  volumetric  flow  rates  correspond  to  ini¬ 
tial  reactant  velocities  of  2.12  mm  s-1  and  6.37  mm  s-1 ,  respectively. 
Two  type  J  and  two  type  K  thermocouples  were  placed  in  10  mm 
increments  down  the  length  of  the  MPB  reactor  channel.  A  National 
Instruments  USB  9211  digital  to  analog  converter  connected  to 
9162  USB  data  acquisition  device  transmitted  the  temperature 
data  to  the  computer  for  temperature  data  collection  and  stor¬ 
age.  The  effluent  multiphase  product  stream  consisting  primarily 
of  liquid  H202  and  H20  and  gaseous  02  flowed  from  the  MPB 
reactor  to  the  liquid  collection  vial.  The  02  gas  proceeded  to  an 
Aalborg  oxygen  flow  meter.  A  National  Instruments  USB  6008 
digital-to-analog  converter  transmitted  the  flow  data  to  the  lap¬ 
top  for  collection  and  storage.  The  program  Labview  Signal  Express 
controlled  both  the  temperature  and  flow  data  recording  pro¬ 
cesses.  At  the  end  of  each  experimental  run,  standard  potassium 
permanganate  (KMn04)  titrations  were  conducted  with  the  liquid 
solution  collected  from  the  MPB  outlet  in  the  sample  vial  to  deter¬ 
mine  the  steady-state  outlet  concentration  of  H202  and  conversion 
[22]. 

MPB  reactor  preparation.  The  results  of  the  MPB  modeling  efforts 
were  used  as  the  basis  of  the  MPB  design  completed  in  a  CAD  pro¬ 
gram.  The  experimental  MPB  reactor  was  fabricated  at  a  precision 
machine  shop  using  type  316  stainless  steel.  Glass  wool  was  plugged 
at  the  inlet  and  outlet  to  keep  catalyst  within  the  reactor  microchan¬ 
nel.  Mn02  catalyst  (Mallinckrodt  Baker:  product  codes  5308, 8329) 
with  particles  ranging  in  diameter  from  60  to  110  |xm  and  an  aver¬ 
age  particle  diameter  of  75  p,m  was  applied  to  the  MPB  reaction 
channel.  The  reaction  rate  constant  per  catalyst  mass  was  estimated 
during  kinetic  experiments  with  the  Mn02  catalyst  particles  and 
provided  the  values  of  the  activation  energy  and  frequency  factor 
used  in  the  FEM  model.  The  differential  reactor  mass  was  recorded 
and  gave  the  Mn02  catalyst  mass,  which  was  0.014  ±  0.002  g  for  all 
experiments. 

3.  Results  and  discussion 

3.1.  Simulated  MPB 

Fixing  the  height  of  the  simulated  reactor  channel  to  0.5  mm 
and  varying  its  length,  a  length  of  5  cm  yielded  42%  conversion  at 
the  outlet  using  50%  (w/w)  H202  reactant  and  an  initial  volumet¬ 
ric  flow  rate  of  0.15  ml  min-1  (initial  fluid  velocity  of  6.37  mm  s-1 ). 
However,  a  temperature  rise  of  92  K  was  simulated  across  the 
microchannel  reactor  lacking  the  extended  surface  area  as  shown 
in  Fig.  3.  Despite  microscale  geometry  (S/V  of  2.3  x  103m2m-3) 
the  simulated  reactor  channel  without  extended  surface  area  failed 
to  prevent  significant  heat  rise.  An  outlet  temperature  exceeding 
373  K  would  lead  to  unwanted  water  vaporization  in  the  effluent 
stream  and  subsequently  increased  pressure,  both  of  which  were 
outcomes  beyond  the  scope  of  the  simulation  assumptions.  This 
finding  reinforced  the  inherent  risk  of  thermal  runaway  and  other 
possible  thermal  effects  due  to  the  exothermic  nature  of  the  H202 
decomposition  reaction. 
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Consequently  the  surface  area  surrounding  the  simulated  reac¬ 
tor  channel  was  extended  to  improve  passive  cooling  of  the 
simulated  MPB  reactor  by  increasing  the  S/V.  The  final  simulated 
MPB  reactor  geometry  included  the  stainless  steel  block  of  the 
dimensions  displayed  in  Fig.  1 A  and  had  a  S/V  of  188  x  103  m2  m-3, 
approximately  two  orders  of  magnitude  larger  than  the  original 
microchannel  S/V.  Fig.  3  also  shows  the  resulting  steady-state  tem¬ 
perature  (I<)  and  H202  concentration  (molnrr3)  profiles  down  the 
length  of  the  simulated  reactor  channel  with  an  extended  surface 
area.  The  temperature  rise  down  the  length  of  the  simulated  reac¬ 
tor  channel  surrounded  by  the  extended  surface  area  was  only  7  K 
with  an  outlet  temperature  of  299  K.  The  extended  surface  area 
significantly  enhanced  thermal  control  of  the  simulated  reaction 
zone  and  ensured  the  outlet  temperature  remained  well  below  the 
boiling  point  of  water  (373  K)  while  achieving  a  conversion  of  38% 
calculated  0.25  mm  from  the  reactor  wall. 

The  outlet  temperature  variation  in  the  radial  direction  ranged 
from  an  average  of  299.53  K  to  a  maximum  of  299.75  K,  and  was 
insignificant  compared  to  the  change  in  the  axial  direction.  The 
outlet  H202  concentration  in  the  radial  direction  ranged  from  an 
average  of  10,852  mol m-3  to  a  minimum  of  6283  mol m-3.  The 
concentration  variation  over  the  outlet  resulted  from  decreasing 


Length  down  simulated  reactor  channel  /cm 

Fig.  3.  Temperature  in  K  and  H202  concentration  in  molirr3  down  the  length  of 
the  simulated  reactor  channel  using  an  initial  reactant  concentration  of  50%  H202 
and  volumetric  flow  rate  of  0.15  ml  min-1  simulated  with  and  without  an  extended 
surface  area  surrounding  the  reactor  channel. 


velocity  of  the  fluid  towards  the  rector  wall  and  subsequently 
greater  residence  times.  There  is  slightly  more  conversion  of  hydro¬ 
gen  peroxide  at  the  outlet  of  the  simulated  MPB  reactor  without 
extended  surface  area  because  of  the  additional  thermal  energy 
available  to  promote  the  reaction. 

Simulations  were  run  at  lower  initial  volumetric  flow  rates 
using  the  MPB  geometry  including  the  extended  surface  area  for 
improved  heat  exchange.  Initial  volumetric  flow  rates  lower  than 
0.15  ml  min-1  (initial  fluid  velocity  of  6.37mms_1)  also  exhibited 
thermal  control  over  the  simulated  reaction  zone  by  limiting  the 
release  of  reaction  heat.  Fig.  4  shows  the  simulation  results  of  the 
three-dimensional  finite  element  MPB  model  using  an  initial  50% 
H202  reactant  and  volumetric  flow  rate  of  0.05  ml  min-1  (initial 
fluid  velocity  of  2.12  mm  s-1 ),  the  minimum  flow  rate  used  in  this 
study. 

A  corresponding  temperature  rise  over  the  MPB  reactor  of  5  K 
and  a  conversion  of  63%  calculated  0.25  mm  from  the  reactor  wall 
were  obtained.  The  outlet  temperature  variation  in  the  radial  direc¬ 
tion  ranged  from  an  average  of  297.72 1<  to  a  maximum  of  297.90  K. 
The  outlet  H202  concentration  in  the  radial  direction  ranged  from 
an  average  of  6505  mol  m-3  to  a  minimum  of  3371  molnrr3.  The 
lower  reactant  volumetric  flow  rate  effectively  increased  residence 
time  in  the  reactor  channel  in  the  axial  direction  causing  greater 
conversion  per  unit  volume  of  reactant.  In  addition  to  increased 
conversion,  there  was  also  decreased  temperature  rise  due  to  the 
lower  processing  rate.  For  the  simulated  conversion  of  63%  at  the 
minimum  0.05  ml  min-1  volumetric  flow  rate  the  corresponding 
oxygen  flow  rate  was  5.8  ml  min-1  at  the  outlet.  Despite  the  lower 
simulated  conversion  of  37%  at  the  maximum  0.15  ml  min-1  vol¬ 
umetric  flow  rate,  the  corresponding  oxygen  flow  rate  would  be 
10.1  ml  min-1  at  the  outlet.  Although  the  smaller  volumetric  flow 
rate  resulted  in  higher  conversion  in  conjunction  with  decreased 
temperature  rise,  the  amount  of  total  oxygen  generated  per  unit 
time  exiting  the  reactor  is  reduced  at  the  lower  reactant  volu¬ 
metric  flow  rate.  Thus  the  simulations,  which  are  based  on  source 
terms  for  the  kinetic  reaction  expressions,  presented  an  inherent 
trade-off  between  conversion  efficiency  and  outlet  oxygen  flow  rate 
dependent  upon  the  initial  reactant  volumetric  flow  rate.  A  future 
parametric  study  of  the  MPB  reactor  is  suggested  to  discover  the 
optimum  initial  volumetric  flow  rate  based  on  the  oxygen  flow  rate 
requirement  for  the  UUV  fuel  cell  system. 
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Fig.  4.  Three-dimensional  temperature  distribution  in  K  (top)  and  H202  concen¬ 
tration  in  molm-3  (bottom)  of  the  simulated  MPB  run  using  an  initial  reactant 
concentration  of  50%  H202  and  volumetric  flow  rate  of  0.05  ml  min-1  and  convection 
coefficient  of  10Wm-2  I<-1. 


3.2.  Experimental  vs.  simulated  temperature 

Using  the  temperature  data  collected  from  the  experimen¬ 
tal  MPB,  the  experimental  and  simulated  temperatures  were 
compared.  Fig.  5  gives  the  average  experimental  steady-state  tem¬ 
peratures  recorded  at  10  mm  increments  down  the  length  of  the 
MPB  reactor  in  comparison  to  the  simulated  steady-state  tempera¬ 
ture  profiles  given  by  the  trend  lines  down  the  center  length  of  the 
simulated  MPB  reactor  using  an  initial  50%  (w/w)  H202  reactant  at 
the  minimum  volumetric  flow  rate  of  0.05  ml  min-1  (2.12  mm  s-1 ) 
and  maximum  volumetric  flow  rate  of  0.15  ml  min-1  (6.37  mm  s-1 ). 


The  temperature  rise  was  less  than  7  K  for  both  the  simulated  and 
experimental  MPB  confirming  improved  thermal  control  upon  inte¬ 
gration  of  the  extended  surface  area. 

Fig.  5  also  shows  that  the  experimental  MPB  readings  were 
within  a  2I<  range  of  those  simulated.  Although  the  simulated 
and  experimental  MPB  temperatures  were  similar,  the  increase  in 
temperature  approaching  the  MPB  outlet  at  5  cm  seen  during  simu¬ 
lation  was  not  replicated  during  experimentation.  This  discrepancy 
was  partially  due  to  uneven  catalyst  loading  in  the  experimen¬ 
tal  MPB  reactor  which  was  exacerbated  during  the  progression 
of  the  multiphase  reaction.  As  the  reaction  progressed,  the  ini¬ 
tially  packed  catalyst  particles  were  dislodged  and  aggregated  into 
irregularly  distributed  catalyst  regions.  The  experimental  temper¬ 
atures  were  lower  than  expected  3  cm  down  the  reactor  channel 
because  a  void  in  the  catalyst  packing  was  observed  at  that  loca¬ 
tion.  Two  and  four  centimeter  down  the  length  of  the  MPB,  larger 
clumps  of  catalyst  aggregated  and  temperature  recordings  were 
correspondingly  higher  due  to  greater  localized  conversion.  Fligher 
temperatures  were  observed  for  both  the  simulated  and  experi¬ 
mental  maximum  volumetric  flow  rate  as  well.  This  was  expected 
because  there  was  more  reactant  available  for  conversion  and  cor¬ 
respondingly  exothermic  heat  generation  at  the  higher  volumetric 
flow  rate.  The  experimental  temperatures  were  consistently  lower 
than  those  simulated  at  both  volumetric  flow  rates.  The  decrease  in 
experimental  MPB  temperatures  relative  to  the  simulation  under 
equivalent  conditions  was  due  to  the  lower  experimental  conver¬ 
sions. 

3.3.  Experimental  vs.  simulated  conversion 

Fig.  6  shows  the  experimental  steady-state  outlet  conversions 
were  lower  than  those  simulated  under  equivalent  conditions. 
Since  the  reactant  dispensed  at  0.05  ml  min-1  experienced  greater 
residence  time,  it  produced  higher  simulated  (63%)  and  exper¬ 
imental  (18%)  conversions  relative  to  the  simulated  (37%)  and 
experimental  (5%)  conversions  resulting  from  reactant  dispensed 
at  0.15  ml  min-1.  Lower  experimental  conversions  resulted  in  the 
low  effluent  02  flow  rates  measured  from  the  experimental  MPB. 
When  the  syringe  pump  dispensed  the  50%  H202  reactant  at 
0.05  ml  min-1 ,  the  average  experimentally  measured  effluent  02 
flow  rate  was  0.84  ml  min-1.  When  the  syringe  pump  dispensed 
the  50%  IT202  reactant  at  0.15  ml  min-1 ,  the  average  experimentally 
measured  effluent  02  flow  rate  was  1.15  ml  min-1 . 

3.4.  Factors  reducing  experimental  conversion 

In  addition  to  unbalanced  catalyst  loading  and  an  inherent 
propensity  for  particle  aggregation  during  reaction  in  the  experi¬ 
mental  MPB  reactor  mentioned  in  Section  3.2,  the  development  of 
02  gas  bubbles  at  active  catalyst  sites  likely  reduced  experimental 


Fig.  5.  Experimentally  recorded  and  simulated  steady-state  MPB  temperatures 
down  the  reactor  length  using  initial  50%  (w/w)  H202  reactant  with  volumetric  flow 
rates  of  0.05  ml  min-1  and  0.15  ml  min-1 . 


Fig.  6.  Experimental  and  simulated  steady-state  outlet  conversions  using  initial 
50%  (w/w)  H202  reactant  and  volumetric  flow  rates  (VFR)  of  0.05  ml  min-1  and 
0.15  ml  min-1. 
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conversion  and  overall  experimental  oxygen  production.  The  obser¬ 
vation  of  O2  bubbles  intermittently  affixed  to  catalyst  particles  in 
the  MPB  microchannel  supports  the  notion  that  active  site  block¬ 
ing  by  02  gas  bubbles  caused  lower  conversion  than  expected.  This 
mass  transfer  limitation  effectively  restricts  continuous  access  of 
H202  reactant  to  active  catalyst  sites  in  turn  decreasing  conversion 
and  subsequently  oxygen  production. 

Observations  of  gas  bubbles  indicated  deviation  from  the 
assumption  of  a  homogenous,  well-mixed  effluent  used  during 
simulation.  The  presence  of  such  gas  slugs  reduces  contact  time 
between  the  reactive  stream  and  catalyst.  This  reduced  contact 
time  essentially  results  in  a  lower  residence  time  and  is  likely  to 
have  contributed  to  the  reduction  in  experimental  conversion.  This 
is  indicative  of  one  of  the  fundamental  challenges  of  multiphase 
reactions  taking  place  in  microscale  reactors. 

Another  factor  influencing  and  likely  suppressing  experimental 
conversion  and  oxygen  generation  in  comparison  to  those  sim¬ 
ulated  under  equivalent  conditions  is  the  type  and  quantity  of 
stabilizer  present  in  the  50%  stock  H202  reactant  used.  Fully  char¬ 
acterizing  the  effect  of  H202  stabilizers  is  challenging  because 
stabilizer  information  remains  proprietary  and  varies  among  sup¬ 
pliers.  Preliminary  experimental  runs  with  a  stock  of  stabilized 
H202  resulted  in  decreased  oxygen  production  over  the  course  of 
an  hour.  Fig.  7  displays  the  experimental  effluent  oxygen  flow  rate 
measured  from  the  MPB  outlet  versus  time  in  minutes  using  the 
stabilized  reactant.  The  decline  in  oxygen  production  is  clear  based 
on  the  trend  line  illustrated.  Although  stabilizers  were  present,  the 
model  is  still  informative  especially  in  regard  to  understanding  the 
thermal  performance  of  the  system  which  is  the  primary  objective 
of  the  study.  To  more  closely  predict  the  conversion,  the  reaction 
rate  law  can  be  modified  to  capture  the  effect  of  the  stabilizer  [23], 
or  simulation  and  experiments  can  be  performed  with  alternate 
noble  metal  catalysts  which  are  less  susceptible  to  degradation 
from  the  stabilizers  [24].  The  effects  of  such  proprietary  stabilizers 
on  decomposition  catalysts  must  be  more  fully  examined  in  future 
studies. 

In  addition  to  improved  characterization  of  peroxide  stabiliz¬ 
ers,  improved  catalyst  application  techniques  would  be  helpful 
in  reducing  the  differences  between  simulated  and  experimental 
results  as  well  as  in  enhancing  overall  oxygen  production  of  the 
microreactor.  For  example,  a  thin  layer  that  acts  as  an  adherent  for 
catalyst  particles  could  be  employed.  Fig.  8  depicts  the  microscale 
reactor  with  a  thin  layer  of  adherent  lining  which  keeps  Mn02  cat¬ 
alyst  particles  fixed  to  the  inside  wall  of  the  channel.  Using  the 
adherent  lining,  which  is  a  carbon  ink  compound  the  decomposi- 


Fig.  7.  Decrease  in  experimental  effluent  oxygen  flow  rate  in  ml  min-1  over 
time  using  stabilized  50%  H202  reactant  with  an  initial  volumetric  flow  rate  of 
0.15  ml  min-1. 


Fig.  8.  Section  of  microscale  reactor  in  which  the  black  adherent  lining,  a  carbon  ink 
compound,  is  used  to  hold  catalyst  particles  to  the  inside  wall  of  the  channel.  The 
resulting  effluent  stream  displays  a  well-mixed  dispersed  somewhat  frothy  phase 
in  which  the  02  gas  bubbles  more  readily  mix  with  the  liquid  solution. 


tion  catalyst  stays  fixed  along  the  wall  of  the  microscale  reaction 
channel.  This  approach  not  only  deters  clog  formation  but  also 
seems  to  promote  better  dispersion  of  the  02  gas  in  the  liquid 
stream  as  the  effluent  appears  more  homogenous  when  the  cat¬ 
alyst  lines  the  wall  of  the  reaction  channel.  The  observation  of  a 
more  readily  dispersed  fluid  in  the  MPB  may  overcome  some  of  the 
identified  drawbacks  to  the  packed  catalyst  approach,  such  as  stub¬ 
born  02  bubbles  immobilized  at  active  catalyst  sites  and  blockages. 
The  exploration  of  optimal  catalyst  application  also  remains  a  topic 
for  future  investigations. 

4.  Conclusions 

With  the  advent  of  microchemical  technology,  there  is  a  tempta¬ 
tion  to  assume  that  even  highly  exothermic  reactions  can  be  made 
to  proceed  nearly  isothermally  simply  by  implementing  microscale 
geometry  [25].  However,  the  significant  temperature  rise  in  the 
simulated  microchannel  reactor  without  extended  surface  area 
emphasized  that  highly  exothermic  reactions  are  not  necessarily 
near-isothermal  even  at  sub-mm  diameters.  Increasing  the  S/V ulti¬ 
mately  achieved  thermal  control  over  the  reaction  zone.  High  area 
heat  exchange  surfaces,  like  the  one  added  to  the  MPB  reactor,  offer 
an  easily  integrated,  passive  mechanism  to  mitigate  the  risk  of  ther¬ 
mal  runaway  in  microscale  reactors.  Although  simulations  were 
shortcoming  in  capturing  the  behaviors  exhibited  in  the  experi¬ 
mental  effluent,  they  indicated  conversions  greater  than  30%  would 
be  possible  under  improved  catalyst  utilization. 

To  increase  experimental  oxygen  production  with  the  MPB, 
a  continuous,  even  distribution  of  catalyst  particles  in  the  reac¬ 
tor  microchannel  is  needed.  Catalyst  application  along  the  reactor 
wall  using  an  adherent  thin  layer  is  one  alternative  that  deserves 
future  investigation.  This  technique  could  offer  more  open  vol¬ 
ume  in  the  microchannel  and  permit  oxygen  gas  to  be  efficiently 
removed,  rather  than  being  trapped  in  the  voids  of  the  packed 
reactor,  potentially  overcoming  the  experimentally  observed  mass 
transfer  limitations  from  the  fixation  of  oxygen  bubbles  at  active 
catalyst  sites. 

Since  commercial  H202  stabilizers  likely  contributed  to  catalyst 
deactivation  of  the  Mn02  catalyst  in  the  work  here,  better  charac¬ 
terization  of  stabilizers  and  their  impact  on  the  reactor’s  catalyst 
requires  additional  investigation.  Although  microscale  reactor  heat 
dissipation  is  easily  improved  passively,  the  MPB  approach  is  sub¬ 
ject  to  inherent  mass  transfer  challenges  that  should  be  addressed 
via  optimization  of  catalyst  application  and  reactor  geometry. 
Resolving  these  issues  will  facilitate  the  realization  of  a  multi-unit 
microchemical  system  for  greater  total  oxygen  production  levels. 
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These  findings  present  a  foundation  for  thermally  controlled  H2O2 
decomposition  microscale  reactors  and  suggest  mechanisms  for 
further  reactor  optimization  to  achieve  increased  oxygen  genera¬ 
tion  rates  for  UUV  SOFC  power  systems. 
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Appendix  A.  Simulation  constants 


Equation  symbol 

Values 

Description 

SI  units 

v0 

2x  10- 

11  to  2  x  10-8 

Initial  volumetric  flow  rate 

m3  s_1 

Vxo 

5x  10- 

5  to  5  x  10-2 

Initial  superficial  fluid  velocity  in  x-direction 

ms-1 

To 

293 

Initial  temperature 

K 

Ea 

20,400 

Reaction  rate  law  activation  energy 

J  mol-1 

A 

2,910 

Reaction  rate  law  frequency  factor 

mol  m-3  s-1  gcat-1 

Ph2o2 

1097 

Average  H202  density  of  solution  in  reactor 

kg  m-3 

Ph2o 

998 

Density  of  water 

kg  m-3 

po2 

1.3 

Average  02  density  between  273  and  350  K 

kg  m-3 

I<e  ff 

5.65 

Effective  thermal  conductivity  of  microchannel 

Wm-1  K"1 

^ambient 

290 

Ambient  temperature 

K 

hconv 

10 

Estimated  natural  convective  heat  flux  for  air 

Wm-2  K"1 

mwp 

0.034 

Molecular  weight  of  H2O2 

kg  mol-1 

mwH2o 

0.018 

Molecular  weight  of  H20 

kg  mol-1 

A  Hrxn 

98,200 

Heat  of  reaction 

J  mol-1 

R 

8.314 

Ideal  gas  constant 

J  mol-1  K1 

cp 

3731 

Average  H202  reactor  solution  heat  capacity 

J  kg-1  K"1 

Mcat 

0.015-0.030 

Mass  of  catalyst 

g 

^effpH20 

7.85  x  10-10 

Average  effective  diffusivity  of  H2O2  into  water 

m2  s_1 

^effpp 

2.71  x  10“12 

Average  effective  diffusivity  of  H2O2  into  H2O2 

m2  s_1 

As 

8000 

Density  of  stainless  steel  316 

kg  m-3 

Cpss 

500 

Heat  capacity  of  stainless  steel  316 

J  kg-1  K-1 

/Css 

16.3 

Thermal  conductivity  of  stainless  steel  316 

Wm-1  K"1 

cpH20 

75.43 

Heat  capacity  of  water 

J  mol-1  K-1 

Cp02 

29.39 

Heat  capacity  of  oxygen 

J  mol-1  K_1 

P0 

1.01  x  105 

Initial  pressure  at  inlet 

Pa 

H 

1.15  x  10“3 

Average  dynamic  viscosity  of  H2O2  solution 

Pas 

Appendix  B.  Simulation  scalar  expressions 

Equation  symbol 

Expression 

Description 

SI  units 

rH202 

krxn  Ch2  0; 

2  x  MCat 

Rate  law  for  H2O2  consumed 

mol  m  3  s_1 

rH20 

rH202 

Rate  law  for  H20  produced 

mol  m  3  s_1 

r°2 

0.5(rn2o 

) 

Rate  law  for  02  produced 

mol  m  3  s_1 

k\xn 

Ae-Ea/RT 

Arrhenius  reaction  rate  constant 

mol  m-3  s-^cat-1 

AConv 

p 

0 

b 

n  1 
P 

0 

P 

0 

b 

Conversion 

1 

Co2 

0.5Xconv 

T0 IT  O2  concentration 

mol  m-3 

'-H202,0 

(l+0.5.XH2o2)OXconv) 

Ch2o 

Ch2o2,o((^h2o,oAh2o2,o) 

^conv )  H20  concentration 

mol  m-3 

*h2o2 

Ch2o2/(Ch2o2  +  Co2  +Ch2o) 

Mole  fraction  H202 

1 

*h2o 

Ch2o/(Ch2o2  +  Q>2  +Ch2o) 

Mole  fraction  H20 

1 

yo2 

Co2/(Ch2o2  +  Cq2  +Ch2o) 

Mole  fraction  02 

1 

sh 

cph2o(373  -  70)rH2 0  +  cpo2 (373  -  T0)r02  Sensible  heat 

Wm-3 

A//Vap 

-49.0047+58754 

Heat  of  water  vaporization 

J  mol-1 

^H202eff 

^effpp(^H202  )  +  ArffpH2o(*H2C 

i )  Diffusion  of  H2O2  into  solution 

m2  s-1 

Appendix  C.  Boundary  conditions 

Boundary 

Boundary  condition  Equation 

Mass  balance 

Inlet 

Initial  H2O2  concentration  C  =  Ch2o2,o 

Outlet 

Convective  flux;  no  diffusive  flux  n(-DeffVCH2o2  =  0)a 

Reactor  walls 

Symmetry 

n(-DeffVCn2o2  +  C\\2o2Vx  =  0) 

Energy  balance 

Inlet 

Initial  temperature 

7=  70 

Outlet 

Convective  flux;  no  conductive  flux  n  ■  (-/CeffVT  =  0) 

Internal  surfaces 

Continuity 

n  (<7i -q2)  =  0]  q{  =  -XeffiV7i  +  p,cPi (7,7,- 

External  surfaces 

Heat  flux 

It  ■  (— keffV7  +  pCpUT )  =  hConv(P  —  ^ambient) 

Momentum  balance 

Inlet 

Inflow  velocity 

Ef  —  vXoi 

Outlet 

Normal  flow/pressure  P  =  Patm 

Reactor  walls 

No  slip 

U=  0 

a  n  represents  the  normal  outward  vector. 
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